Incoherent matter-wave solitons: Mutual self-trapping of 
a Bose-Einstein condensate and its surrounding thermal cloud 
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We show that a Bose-Einstein condensate and a portion of its surrounding thermal cloud can ex- 
hibit mutual self-trapping, supported by the attractive particle interactions, and not by the external 
confinement. This type of dynamics is characteristic of composite random-phase solitons. 



The physics of quantum-degenerate, interacting bose 
gases closely resembles the behavior of light in non- 
linear media. The dynamics of Bose-Einstein conden- 
sate (BEC) at zero-temperature is within the Gross- 
Pitaevskii (GP) mean-field theory described by the non- 
linear Schrodinger equation (NLSE) for the condensate 
order parameter. The same equation describes the evo- 
lution of coherent light in nonlinear Kerr medium. This 
analogy has opened the way for the field of nonlinear 
atom optics |2j with striking demonstrations of famil- 
iar nonlinear optics phenomena such as four wave mix- 
ing , superradiant Rayleigh scattering Q , and matter- 
wave amplification 0, carried out with matter- waves. 
One suchphenomenon is the formation of matter-wave 
solitons H H [H Elljrljri 0i0. Experi- 
mentally, dark solitons [3, Ui| an( l bright gap solitons 
|l6j were observed in BECs with repulsive interactions, 
whereas bright solitons fl^fl4| were demonstrated in sys- 
tems with attractive interactions. These experimental 
results are augmented by extensive theoretical work in- 



and dark ^lj matter- 



eluding predictions on bright 

wave solitons, lattice solitons [12j , and soliton trains [If 
To the best of our knowledge, all previous theoretical 
efforts on matter-wave solitons have utilized the zero- 
temperature GP mean-field theory. However, in a realis- 
tic system, elementary excitations arising from thermal 
and/or quantum fluctuations are always present, and the 
BEC dynamics may be considerably affected by the mo- 
tion of excited atoms around it (thermal cloud), giving 
rise to new nonlinear matter-wave phenomena. 

Here we present an example of such novel phenomena, 
and show that a BEC and a portion of its surrounding 
thermal cloud can exhibit mutually self-trapped motion. 
This motion is achieved via attractive interactions be- 
tween particles, and not by the external confinement. We 
emphasize that the finite-temperature self-trapping pro- 
duces a truly novel type of matter-wave solitons, where 
localization is attained not only in spatial density but 
also in spatial correlations. These self-trapped incoherent 
matter-waves are analogous to c omp osite random-phase 
(incoherent) optical solitons |ljl Ilia, LUj ; one component 
is the BEC while the other is a part of its thermal cloud. 
An important result of this work is that the established 



analogy between zero-temperature BECs and coherent 
nonlinear optics can be elevated to the analogy of inco- 
herent light behavior in nonlinear media and BECs at 
finite-temperatures. 

The incoherent matter-wave self-trapping is illustrated 
within the Hartree-Fock (HF) approximation. First, we 
solve the static HF problem [2(1 EH E3 | describing Bose 
gas with attractive interactions in a confining harmonic 
potential. We consider the range of (higher) tempera- 
tures where excitations are mostly particle-like, and HF 
solutions are approximately equal to the solutions ob- 
tained using the Hartree-Fock-Bogoliubov (HFB) theory 
in the Popov approximation [il]]. We focus on solutions 
with significant population of the first-excited state. In 
order to establish mutual-sclf-trapping, we turn-off the 
external harmonic confinement, and evolve the BEC and 
the thermal cloud within the time-dependent HF the- 
ory (TDHF) jH Ullill. The signature of mutual self- 
trapping is a slow separation of the total density profile 
into two humps, which propagate almost in parallel, and 
are then pulled back to almost recover the initial den- 
sity profile. This self-trapped motion is compared to the 
evolution of precisely the same initial state when both 
the trap and the interparticlc interactions are turned off 
simultaneously, exhibiting fast matter-wave dispersion. 

We consider a system of N interacting bosons placed 
in a quasi one-dimensional (Q1D) cigar-shaped harmonic 
potential V ex t(x, y, z) = (u> x x 2 + u>±y 2 + u>j_z 2 )/2, where 
u>± ^> lo x denote the transverse and the longitudinal fre- 
quencies of the trap, respectively. The interparticlc in- 
teraction is approximated by the Q1D contact potential 
V(xi - x 2 ) = 9idS{xi - X2), where g 1D = -2h 2 /ma 1D , 
a\D ~ — a^/a^D is the effective ID scattering length 
[2^, [23, m is the particle mass, a± = ^JKJmujZ is the 
size of the lowest transverse mode, while a^o is the 3D 
scattering length. At finite temperatures, this system 
can be describe d by the Hartree-Fock Bogoliubov (HFB) 
pairing theory [21]. We are interested in the regime of 
higher temperatures, where the thermal cloud is suffi- 
ciently large; in this limit, anomalous pairings are neg- 
ligible, the excitations are mostly particle-like, and the 
HFB formalism can be approximated with the HF the- 
ory [iljl . In our case, the system is initially at thermal 
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equilibrium. The density of the condensate \(ip(x,t))\ 2 , 
and the noncondensed particles W (x, i)il)(x , t)) are cal- 
culated within the static HF approximation |20ll2lll2^ |: 

H SP 4 S) + g 1D [n{ s \x) + 2n[ s \x)]4 s) = ec4 s) (x), (1) 



H sp uf + g 1D {2ni s \x) + 2n[ s \x)}u < j s) = e jU ^(x). (2) 

Here, H sp = -^gjT + \mLo 2 x 2 \ (f>o\x) is the conden- 
sate wavefunction, with eigenvalue eo; u^ s \x) are the 
wavefunctions of the excited states, with eigenvalues tj, 

j = 1,2,...; nc (x) = iV c |(/f)Q S ^(x)| 2 is the condensate 
density, where N c is the number of condensed particles; 

J2j Aj|u^| 2 is the density of the excited par- 



0)/ i 



tides; the number of particles populating the jth ex- 
cited state, Nj, is determined by the Bose distribution, 
Nj = [exp( c £~t ) — 1] — 1 - The chemical potential p is 
set by the constraint N = N c + J2j Nj ; all wavefunc- 
tions are normalized to unity. We solve Eqs. and @ 
sclf-consistcntly, and use the static solution as the initial 
condition to study dynamics without confinement. 

The dynamics is studied within the TDHF approxima- 
tion, which involves the coupled evolution of the conden- 
sate wavefunction <f>o(x,t) = (ip(x,t)}, and the correla- 
tion function p(x x , x 2 , t) = (^(x 2 ,t)i)(x 1 ,t)) [HlilE^: 



ih 



dp 



dcj)o(x,t) 



Of 



H sp 4>o + giD[n c (x,t) + 2nt(x,t)]<j>o, (3) 



ih^ = [H sp (l) - H ap (2)]p + 2g 1D [n c (xi,t)+n t (xi,t) 



n c (x 2 ,t) - n t {x 2 ,t))p(xi,x 2 ,t), 



(4) 



where n c (x,t) = N c \4>q(x, t)\ 2 and rit(x,t) = p(x,x,t). 
The correlation function p is a part of the single particle 
density matrix (ip^ (x 2 ,t)ij}(xi,t)) = <t>o(x 2i t)<f>o(xi,t) + 
p(xi,x 2 ,t) corresponding to excitations. Equations of 
motion © and (0} for the condensate wavefunction and 
the excitations describe the evolution of the BEC and 
the thermal cloud even outside of equilibrium Ini- 
tial conditions at t = are <fio(x,t = 0) = <p Q (x) and 
p(xi,x 2 ,t = 0) = Y^jNjU*^(x2)Uj(xi), and the evo- 
lution is performed without the external potential. How- 
ever, rather than to use Eqs. (0 and (0J, we shall use 
a fully equivalent, but numerically more convenient ap- 
proach. The solution for p(xx,x 2 ,t) may be constructed 
from p(x\, x 2 , t) = NjUj(x2,t)iij(xi,t), where func- 
tions Uj (x, t) evolve according to an infinite set of coupled 
equations, 



duj(x, t) 
dt 



H sp Uj + g 1D 2[n c (x,t) + n t (x,t)]uj, (5) 
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FIG. 1: (color online) Solution of the static Hartree-Fock 
equations, (a) The total density (solid line), the bell-shaped 
condensate wavefunction (dotted line), and the dipole-shaped 
first excited state (dashed line), (b) The (mean-field) poten- 
tial seen by the condensate (solid line) , by the non-condensed 
particles (dotted line), and the external harmonic potential 
without the mean field (dashed line). Horizontal lines from 
the lowest one up depict the chemical potential p, the con- 
densate eigenvalue eo, eigenvalue ei of the dipole-mode, and 
eigenvalues e2, and e$, expressed in units of hu>±. 



f (x,i)| : 



Note that when 



where n t (x,t) 
the time-dependence of the wavefunctions is sim- 
ply 4>o(x,t) = (a;) exp(— ieot/K) and Uj(x, t) = 
u^ s \x) exp{—iejt/K), the equations of motion JSJl and (JHJ 
reduce to the static HF equations JIJ and 

In what follows we present results of a numerical calcu- 
lation based on the described HF formalism, demonstrat- 
ing the mutually self-trapped motion of the BEC and a 
portion of its thermal cloud. The parameters of the cal- 
culation are chosen to resemble the experimental param- 
eters of Ref. 0. We consider TV = 2.2 10 4 7 Li atoms in 
a harmonic trap with uj± = 4907 Hz (aj_ = \Jh/muJi_ ~ 
1.35 jLtm), and u> x = 35 Hz (a x = y h/moj x 



» 16.0 pm). 

The 3D scattering length cl^d = —3.1 10~ n m corre- 
sponds to a nonlinear parameter of iV|a3o| « 0.68 urn , 
and is tunable by the Feshbach resonance technique [l3j . 
The temperature is fc^T/liWi = 16. The notable differ- 
ences from the experiment of Ref. 0| are in the longi- 
tudinal trapping frequency, which is an order of magni- 
tude smaller in our calculation, and in the temperature, 
which is here chosen to produce a sufficiently large ther- 
mal cloud and to ensure the validity of the HF formalism. 

While the temperature kgT in our simulation is higher 
than the transverse level spacing tuv \ whereas a 'true' ID 
geometry calls for fc^T < huj± |26J,|2jJ, the use of a Q1D 
formalism is still justified because the first lo x /lo± ~ 140 
states are essentially ID (they are in the lowest state of 
the transverse Hamiltonian). As shown below, only a 
few of the lowest excited states actively participate in 
the self-trapping process. Therefore, a proper inclusion 
of the transverse dimension in the calculation would lead 
to some rescaling of the parameters, but would not influ- 
ence the self-trapping process observed in our quasi-lD 
calculation. Moreover, the simulations as well as the ex- 
periment of |l3j , are far from the Tonks-Girardeau regime 
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FIG. 2: The evolution of the system (a) with interactions 
present, and (b) without interactions. In both cases, ini- 
tial conditions were simply the static HF solution from Fig. 
while external potential was turned off during evolution. 
Graphs show the total density profiles at equally spaced in- 
tervals from t = up to t = 1.2 x 2n/u> x ; the x-coordinate is 
£ = x/a±. 



of impenetrable bosons [26j ; our calculations are all in the 
weak interaction regime N\aiu\/a x ~ 10 7 3> 1, thus jus- 
tifying the use of a mean-field approach (we note paren- 
thetically that the weak interaction regime for a ID gas is 
attained 'counterintuitively' at high densities |26jV The 
stability of the confined, attractively interacting conden- 
sate against collapse 0, was numerically verified by 
evolving it with random initial noise on top of all equi- 
librium modes; the stability is underpinned by the use of 
parameters resembling the experiment |l3j . 

The system is initially in equilibrium. Fig. ^a) illus- 
trates the total density (x) + n\ (x) of the stationary 
HF calculation. The total density profile (solid line) is 

double humped, offering a clear signature of the signifi- 
es') 

cant population of the first excited u\ state, which has 
a dipole-like spatial profile (dashed line). The conden- 
sate fraction is N c /N 0.24, while the population of 
the first excited state is Ni/N m 0.093. The excitation 
energies and the mean-field potentials affecting the con- 
densate and the non-condensed cloud arc plotted in Fig. 
Bib). The eigenvalue e\ of the first excited state is below 
zero, along with the condensate eigenvalue eo and the 
chemical potential fx (p, — ► eo when T — > 0). Thus, a 
large fraction of the thermal cloud [Ni/(N - N c ) ps 0.12] 
is self-localized jointly with the condensate, rather than 
localized by the external harmonic potential. The same 
conclusion is also inferred from the double- well structure 
of the mean-field potentials [Fig. [Hb)]. 

The static evidence for mutual self trapping is con- 
firmed by studying the evolution of the system, once 



the trapping potential along x is suddenly turned off, 
as in the experiment of Ref. 0]. The system is sud- 
denly taken out of equilibrium and consequently starts 
to evolve. We simulate the dynamics by solving Eqs. © 
and JSJ with the standard split-step Fourier technique. 
In the spirit of Ref. [l3T | , we compare the x-unconfined 
dynamics of the system in the presence of interparticle 
interactions [Fig. Ufa)] to its time evolution when both 
the confinement in x and the interactions are turned off 
[Fig. Hfb)]. In the absence of interactions, we clearly ob- 
serve a fast dispersion of the total density. After approx- 
imately td ~ 1/3 x 2tt/lo x , there are hardly any particles 
left within the window shown in Fig. |2[b). In contrast, 
when interactions are present, we observe self-trapped 
motion [Fig. |^a)]. The two humps begin to separate, 
because the trapping potential which provided a balance 
to the kinetic energy term is no longer present. How- 
ever, due to the attractive particle interactions, the two 
humps separate very slowly and move almost in paral- 
lel [Fig. EJa)]. Subsequently, the two humps are pulled 
back, and the initial density roughly recovers its initial 
appearance. We emphasize that a significant portion of 
the atoms within the self-trapped entity (Ni/N c ~ 0.38) 
belong to the thermal cloud, and that they strongly affect 
the BEC dynamics. Atoms from higher excited states 
(j > 10) disperse quickly even with interactions present, 
and are essentially spectators of the self-trapped motion. 

After the confining potential is turned off, the con- 
densate would be completely depleted after some time- 
period, which is proportional to the collision time r c , 
because relaxation is attained through collisions. The 
TDHF formalism is clearly inadequate to depict this de- 
pletion since it conserves all populations including the 
condensate fraction. However, since in the weak interac- 
tion regime the collisional time is much longer than the 
trap period t c uj x 3> 1 |30l |. relaxation is slow compared 
to characteristic trap times. We can therefore safely use 
the TDHF method to carry out simulations lasting up 
to several dispersion times Tdisp ~ 1/3 x 2tt/lu x <C t c in 
order to demonstrate self-trapping. 

The self-trapped entities presented here represent par- 
tially coherent matter waves. Thus, it is important to 
study the complex degree of coherence of these matter- 
waves, p(x,x',t) = p(x, x',t)/y/ p(x, x, t)p(x',x', t). Fig. 
[3]shows the evolution of \p{x, x',t)\ corresponding to our 
self-trapped entity. We see that spatial correlation is fi- 
nite. During the initial stage of the evolution \p(x, x', t)\ 
broadens indicating the increase of coherence, but after 
this initial stage it stays practically unchanged. This cor- 
responds to the fact during the initial stage of the evolu- 
tion, a portion of the thermal cloud that is not mutually 
trapped with the BEC disperses, which initially increases 
the coherence. We emphasize that for zero-temperature 
GPE solitons, the pair correlation function factorizes as 
p(x,x') = <jf{x)(j){x'), which yields p,(x,x') = 1, corre- 
sponding to coherent matter- waves. Our incoherent self- 
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FIG. 3: (color online) The complex degree of coherence 
\/j,(x,x' ,t)\ of a partially-coherent, self-trapped matter-wave 
at times t = (black dashed line), t = QAAt (red solid line), 
t = 0.8At (green dot-dashed line), and t — 1.2 At (blue dotted 
line); £ = x/ax, x' = 0, and At = 2tt/lox. 



trapped matter-waves are thus rather special in that they 
correspond to localization of entropy and spatial correla- 
tion, as well as to localization of density. 

Before closing, we return to the analogy between the 
propagation of incoherent light in nonlinear media, and 
the behavior of BECs at finite temperatures. Incoher- 
ent light can be described by the mutual coherence func- 
tion B{x\,X2,z) = (E*(x 2 , z,t)E(xi, z,t)) 0, where 
E(x,z,i) is the randomly fluctuating field. The evolution 
of B(xi,X2, z) along the propagation axis z is described 
by an equation equivalent in structure to Eq. J2J de- 
scribing the evolution of correlations p(xi,X2,t) in time 
(see Ref. 0). Moreover, the TDHF equations JHJ are 
analogous to the Manakov equations that describe inco- 
herent light in nonlinear media fl8t Il9j . The analogy 
is not complete, due to the fact that the Bose wavefunc- 
tion is affected by a different mean-field than the thermal 
cloud, while in optics it is usual (but not the rule) that 
all fields 'see' the same nonlinear change in the index of 
refraction. Furthermore, in the full HFB approximation, 
the U(l) symmetry is broken, which yields phenomena 
not encountered in incoherent nonlinear optics. Never- 
theless, we expect that many nonlinear phenomena with 
incoherent light will find its counterpart in BECs at finite 
temperatures. 

In conclusion, we have demonstrated that a BEC and a 
portion of its surrounding thermal cloud can exhibit mo- 
tion analogous to composite random-phase (incoherent) 
optical solitons. The relation between BEC and nonlin- 
ear optics is thus elevated to the analogy between nonlin- 



ear incoherent optical waves and nonlinear matter-waves 
at finite-temperature. The predicted incoherent matter- 
wave structures represent novel correlation solitons which 
resemble localized second-sound entropy waves. Work is 
underway to go beyond HF approximation and study the 
role of pairing in lower temperature soliton dynamics. 
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